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Abstract. The celebrated problem of a non- homogeneous sphere rolling over a horizontal 
plane was proved to be integrable and was reduced to quadratures by Chaplygin. Applying 
the formalism of variational integrators (discrete Lagrangian systems) with nonholonomic 
constraints and introducing suitable discrete constraints, we construct a discretization of the 
£S| . n-dimensional generalization of the Chaplygin sphere problem, which preserves the same first 

integrals as the continuous model, except the energy. We then study the discretization of 
the classical 3-dimensional problem for a class of special initial conditions, when an analog 
of the energy integral does exist and the corresponding map is given by an addition law on 
elliptic curves. The existence of the invariant measure in this case is also discussed. 
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1 Introduction 



The Chaplygin problem on a non-homogeneous sphere rolling over a horizontal plane without 
slipping is probably one of the best known integrable systems of the classical nonholonomic 
mechanics. Although being non-Hamiltonian in the whole phase space (see [2]), the equations 
of motion possess an invariant measure, which is a rather strong property putting them rather 
close to Hamiltonian systems. 

It is even more interesting that the Chaplygin sphere appears to be an algebraic integrable 
system in the sense that generic level varieties of the first integrals are open subsets of 2-di- 
mensional Abelian varieties and, after an appropriate time reparameterization, the phase flow 
becomes a straight line uniform on them, [3 [8]. 

Note that a Lax pair with a spectral parameter for the Chaplygin sphere, which could provide 
all the constants of motion, is still unknown and, probably does not exist. Hence, one cannot use 
the powerful method of Baker- Akhieser functions to find theta-function solution of the problem 
or to construct its integrable discretization by applying a Backlund transformation, as described 
in [12\ 17]. or by factorization of matrix polynomials. 

Contents of the paper. We briefly recall the equations of motion of the Chaplygin sphere, 
as well as their n-dimensional generalization. Then we construct a discretization of the prob- 
lem by applying the formalism of variational integrators with nonholonomic constraints recently 
developed in [H Q3]. Apart from the nonholonomic distribution on the tangent bundle TQ of 
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the configuration manifold Q, the formalism requires introducing discrete constraints on Q x Q 
which, in a certain sense, must be consistent with the continuous ones. As one can see in the li- 
terature [H [TOl [T5] , the choice of discrete constraints is crucial for the behavior of the discretized 
nonholonomic system: the latter may inherit the main properties of its continuous counterpart, 
or may not. As an example, we note that although continuous systems with stationary nonholo- 
nomic constraints possess the energy integral, almost all their known discretizations do not enjoy 
this property (see e.g., [5]). Nevertheless, for a class of discretizations considered in [10] there 
exists a natural method to specify discrete constraints which ensures the exact preservation of 
energy. 

Our choice of discrete nonholonomic constraints on E(n) x E(n) that mimics the condition 
of the sphere rolling without slipping over a horizontal plane allows us to construct a map 
discretizating the Chaplygin sphere problem, which has the form of a momentum conservation 
law and therefore preserves all the momenta integrals of the original system. 

We then consider the discretization of the classical 3-dimensional problem for the case when 
the discrete angular momentum is vertical. In this special case the structure of the map is 
reminiscent to that of the Veselov-Moser discretization of the Euler top on SO(3) \19\ I16j and 
an analog of the quadratic energy integral does exist. 

This implies that the invariant manifolds of the discretization map are elliptic curves £ and 
the map is described as an addition law. However, in contrast to most known integrable discrete 
systems, in the discrete Chaplygin sphere the translation on £ depends not only on the constants 
of motion but also on the point on the curve. Thus, in order to find the explicit solution, we 
arrive at a rather difficult problem of reparameterization of £ or its real part, which would 
make the translation constant. We notice that this problem is equivalent to the problem of the 
existence of an invariant measure of the map. 



2 The Chaplygin sphere and its multidimensional generalization 

Following Chaplygin [3], consider a dynamically non-symmetric sphere with inertia tensor J, 
radius p, and mass m rolling without slipping over a horizontal plane. Assume that the mass 
center and the geometric center C of the sphere coincide. Let 7 = (7i,72,73) T be the vertical 
unit vector and u> = (07, W2, W3) T , v 6 M 3 be respectively the angular velocity of the sphere and 
the velocity of its center in the moving frame. The condition of non-slipping of the point P of 
contact of the sphere with the horizontal plane is 

v + pojx-f = 0. (2.1) 

Here and below, x denotes the vector product in M 3 . On the configuration space of the problem, 
the Lie group E(3), these equations define nonholonomic constraints: for any two positions of 
the sphere on the plane there exists a linking trajectory that satisfies (|2.ip . 

Under these constraints the equations of motion can be reduced to the following closed system 
for uj, 7 G2 

KCj = kuo x uj + (A _1 7, (Acj) x 01)7, 7 = 7x0;, (2.2) 
F 

where 

D = mp 2 , A = J + DI, F= (J 7 ,A~ 1 7 ) = 1 - £>(7, A _1 7>, 



1 We use our proper notation for the variables and parameters of the problem, which may differ from that used 
in the original paper [3]. 
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I being the identity 3x3 matrix. Let K = Auj — Z) (0^,7)7 be the vector of the angular 
momentum of the sphere with respect to the contact point P. Then the system (|2.2p also 
admits representation in form 

K = Kxuj, 7 = 7x0;. (2.3) 

Hence, like 7, the momentum K is fixed in space. As a result, the system possesses four 
independent first integrals 

(7,7) = 1, (K,j) = (Ju n ) = h, (K,K)=n, (K,u)=l. (2.4) 
Here the last integral is the kinetic energy of the sphere and, since 

it can be rewritten in terms of K, 7: 

(K,A- 1 ^)-|( 7 ,A- 1 K) 2 . (2.5) 
In addition, the system (|2.2p possesses an invariant measure 



VVtTA^t) du! A • • • A dj 3 , (2.6) 

hence, by the Euler-Jacobi theorem (which is also often refereed to as the Jacobi last multiplier 
theorem, see e.g., |2Uj). it is integrable by quadratures and its generic invariant varieties are 
2-dimensional tori. 

There are two special types of the initial conditions, when the equations of motion are simp- 
lified. In the first case the momentum K is horizontal, (K, 7) = 0. As shown by Chaplygin, after 
a time reparameterization and introducing spheroconic coordinates on S 2 = (7, 7) the variables 
separate and the system reduces to hyperelliptic quadratures. Theta function solution for the 
unreduced Chaplygin system in this, as well as in the generic case, was obtained in [8]. 

The second special case is described below. 

The case of the vertical momentum K. As noticed in [3] and as follows from the energy 
integral in (|2.4p . under the special initial conditions K = /17, h = const one has (K,u) = 
h(j,u)) = I, i.e., 

. I h 2 + Dl 

(uj, 7) = — = const, and Au; = 7. 

h h 



As a result, the first vector equation in (12. 2p takes the form of the Euler top equations 

Au = Alo x u, (2.7) 

and integrals (|2.4p reduce to 

(u, Auj) = I, (Au>, Auj) = k 2 , I, k = const. 

Hence, for almost all initial conditions K = /17, the variables Wj, 7 j are elliptic functions of the 
original time t and the solution of the reduced system is periodic. 

Multidimensional generalization. Now, following [9] , consider the generalized Chaplygin 
problem on an n-dimensional ball rolling without slipping on an (n— l)-dimensional 'horizontal' 
hyperplane H in W 1 . The configuration space for the ball is the Lie group E(n) = (R, r), where 



4 



Yu.N. Fedorov 



R G SO(n), r G IR n are the rotation matrix of the sphere and the position vector of its center C. 
For a trajectory R(t), r(t), define the Lie algebra element (w,vc), where 

uj = R- l R g so(n), v c = R~ l r G M. n (2.8) 

are respectively the angular velocity and the velocity of C in the frame attached to the sphere. 

Let now 7 G M n be the unit vector orthogonal to the hyperplane H and directed 'upwards', 
i.e., from H to the center C, and, as above, p be radius of the ball. Then the condition for the 
sphere rolling without slipping that generalizes (|2.ip reads 

v c + pwj = 0. (2.9) 

One can show that this vector constraint determines a non-integrable distribution on the tangent 
bundle T E(n), which is neither left- no right-invariant with respect to the action of E(n). 

Next, introduce the angular momentum of the sphere in the body frame with respect to the 
center C, 

M = Jlo + wJ G so*(n), (2.10) 

J being a constant diagonal n x n matrix. Then, as shown in [9], the motion of the sphere is 
described by the following Euler-Lagrange equations 

M + [w, M] = pF A 7, m(v c + w v c ) = F, 7 + ^7 = 0, (2.11) 

where F G W 1 is the reaction force acting at the point P of contact of the sphere with the 
hyperplane H. Here all the vectors are considered in the frame attached to the ball. 

Now, differentiating the constraints (|2.9p and using the second equation in (|2.1ip . we find 
F = mpojj. Then (|2.1ip gives 

M + [u,M]=D(cbT + Tu), f+[w,r]=0, 

where D = mp 2 , T = 7 <g> 7 = 77^. Then this system can be represented in the following 
compact commutative form that generalizes (|2,3|) 

k + [w,K] = o, r + [w,r] = o, (2.12) 

where 

K = Ju) + ujJ + D{Tuj + ujT) G so*{m) (2.13) 

can be regarded as the angular momentum of the ball relative to the contact point P. 

It is natural to introduce the nondegenerate inertia operator A.: so(n) — > so(n)* such that 
K = Aw. Since A is nondegenerate, equations (I2.12|) give a closed system for the variables w, 7. 

As follows from the form of (|2.12p . Ad^K is a constant tensor in the space. Hence this system 
has the following set of integrals 

trK s , tr (KT 1 ), tr V = 1, s = 2,4,6,..., I G N. (2.14) 

naturally, the system also possesses the energy integral 

H = --tv(wAw) (2.15) 



and, as shown in [9], the invariant measure pdw A dj with density p = VdetA, which depends 
on the components of 7 only. 

In the classical case n = 3, under the isomorphism between so(3) and ]R 3 , 

Wij = EijkWk, K.;j = EijkKk 

the integrals (|2.14p . (|2,15p and the measure transform to (|2,4p and (|2.6p respectively. 

Various properties of the Chaplygin sphere and its n-dimensional generalization were studied 
in [BJ UT\ [51 [TS], although the integrability or nonintegrability of this generalization was not 
proved. 
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3 Discrete mechanical systems with nonholonomic constraints 

Before describing the discrete setting, let us recall that a continuous nonholonomic Lagrangian 
system is a triple (Q,L,TJ), where Q is a smooth n-dimensional configuration space, L : TQ — > 
R is a smooth function called the Lagrangian, and D C TQ is a /c-dimensional constraint 
distribution. Let q = (q±, . . . ,q n ) be local coordinates on Q. In the induced coordinates (q,q) 
on the tangent bundle TQ we write L(q,q). It is assumed that the map 

BT 

— : TQ -> T*Q 
oq 

is invertible. A curve q(t) G Q is said to satisfy the constraints if q(t) G T> q n\ for all t. The 
equations of motion are given by the Lagrange-d'Alembert principle: 
b 

L(q i ,q i )dt = 0, (3.1) 

where 5q(t) G T) q u\ for t G (a, 6) and 5g(a) = = 0. This principle is supplemented by the 
condition that the curve itself satisfies the constraints. Note that we take the variation before 
imposing the constraints. 

Assuming that the constraint distribution is specified by a set of differential forms A 3 (q), 
j = 1, ...,s < n, 

V = {q£TQ\ {A>(q),q)=0, j = 1, . . . , s}, (3.2) 
the principle (|3.1|) implies the Euler-Lagrange equations with multipliers Xj 

d dL 8L sr~^ > . , ; , ^ , 
-^ttt - = ? ^i AJ (q • (3-3) 
dt dq dq ^ 3 w y 1 

Coupled with (|3.2p . they give a complete description of the dynamics of the system. 
Note that equations f|3.3[) conserve the energy 

The discrete Lagrange d'Alembert principle and equations. Let, as above, Q be 
a smooth manifold. According to [3], a discrete nonholonomic mechanical system on Q is 
defined by three ingredients 

(1) a discrete Lagrangian L: Q x Q -> R; 

(2) an (n — s)-dimensional distribution D on TQ (given by equations (|3.2p ); 

(3) a discrete constraint manifold T>d C Q x Q, which has the same dimension as P and 
satisfies the condition (q, q) G P,i for all q £ Q. 

The dynamics is given by the following discrete Lagrange-d'Alembert principle (see (JJ), 
JV-l 

^ (Z>iL(g fe , g fe+1 ) +D 2 I.(gjb-i,5A i ))% = 0, 5q k eV qk , (q k , q k+1 ) G 2? d . 
fc=o 

Here L>iL and -D2L denote the partial derivatives of the discrete Lagrangian with respect to the 
first and the second inputs, respectively. 

The discrete constraint manifold is specified by the discrete constraint functions 

Fj(q k , qk+i) = 0, j = 1, . . . , s, (3.4) 

which impose the above restriction on the solution sequence {(q k ,q k +i)}- 



6 



Yu.N. Fedorov 



Remark 1. If the discrete Lagrangian L is obtained from a continuous one, L(q,q), via a dis- 
cretization mapping ^: Q x Q — > TQ defined in a neighborhood of the diagonal of Q x Q, i.e., 
L = Lo$, then the variety T>d must be consistent with the continuous distribution T>. Namely, 
Vd is locally defined by the equations A 3 ' o ^ = 0, j = 1, . . . , s. 

In particular, if the distribution T> C TQ is integrable, i.e., it defines an (n — s)-dimensional 
integral submanifold M C Q, then T>d C Q x Q is just the product x A/". 

We emphasize that, in general, the discretization mapping is not unique and hence there 
are many ways to define the discrete Lagrangian L and the discrete constraint manifold T>d for 
a given nonholonomic system (Q,L,2?)o 

The discrete principle implies that the solutions of a discrete nonholonomic system are repre- 
sented by sequences {(q k , q k +i)} that satisfy the discrete Lagrange-d'Alembert equations with 
multipliers 

s 

D 1 L(q k ,q k+1 ) + D 2 h(q k _ 1 ,q k ) = ^\ k j A j {q k ). (3.5) 

The multipliers are determined from the discrete constraints (|3.4p . but, in general, not 
uniquely. Hence, the map Q x Q i— ► Q x Q defined by (|3.4f) . fl3.5f) is multi- valued. 

Remark 2. According to [4], equations (13. 5p introduce a well-defined mapping (q k -i,q k ) i— > 
(gfc, (?fc+i), if the (n + s) x (n + s) matrix 

(D x D 2 L d (q k ,q k + x ) A l (q k ) ■■■ A s (q k )\ 

D2^i(qk,qk+i) o ••• o 

V D 2 F s (q k ,q k+1 ) ••• / 

is invertible for each (q k , q k +i) in a neighborhood of the diagonal of Q x Q. 



4 A discretization of the Chaplygin sphere 

Now we apply the above approach to discretize the generalized Chaplygin sphere problem 
on E(n). The trajectory of such system is a sequence (R k ,r k ), k G Z. We choose the dis- 
crete Lagrangian in form 

L = ^Tr{R k JRl +l ) + j (Ar k , Ar k ), Ar k = r k+l - r k . (4.1) 

Here the rotational part ^Tx(R k JR\" +1 ) coincides with the Lagrangian of the Euler top on SO(n) 
introduced in |19j . 

Note that, in view of (|2.8|) . the continuous constraints (|2.9p can be rewritten as 

r + pRR T ^ = 0, (4.2) 

where 7 is the unit normal vector in the fixed frame (without loss of generality we can set it to 
be (0,...,0,1) T .) 

Then, in view of (|4.ip . (|4.2p . m our case the discrete Lagrange-d'Alembert equations (|3.5p 
take the form 

R k+1 J + R k ~\J = A k R k + pj* fR k , (4.3) 

2 An alternative approach to the discretization of nonholonomic systems based on a modification of canonical 
transformations was proposed in [14] . 
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m(r k+ i - r k ) - m(r k - r k -i) = f , (4.4) 

where A^ is the symmetric matrix Lagrange multiplier and /* = . . . , f^) T , is the vector 
multiplier corresponding to the constraints 



In view of the property -Rfc-RjT = I> equations (|4.3p yield 
i? fc +i Ji# + R^JRl = A fe + p/*^, 

i? fc Ji^ +1 + Ji?Li = A£ + P7/1- (4.5) 



Following [19], we introduce the discrete angular momentum of the sphere with respect to its 
center C in space 

m k = Rk+iJRk ~ RkJRl+i G so(n). 
Then (|4"3j) . (fO]) give rise to 

m fc = m fc _i + p(/V - 7/ feT ), m(Ar* - Ar fe _i) = /*. (4.6) 

Next, like in [19], introduce the discrete momentum of the sphere in the body frame 

M k = Ad Rk m fe = n k J - ml, n k = RlR k+l £ SO(n), 

where Q k is the finite rotation in the body frame. In the continuous limit £l k \— > I+eo;, uj G so(n), 
it transforms to the momentum (|2.10p . 
Then the system (|4.6p gives 



M k = Ci*4-i^-i + p(Rk f k j T Rk - RhkFR 



rk 



m{Ar k - Ar fe _i) = f , 
that is, 

M k = Ci*4-ifin + pF k A 7fc , ™(Ar fc - Ar fc _i) = /* , (4.7) 

where 

F fc = 2#/* lk = Rli (4.8) 

represent the vectors 7 in the body frame. Equations (|4.7p can be regarded as a discrete 
version of the equations of motion of the Chaplygin sphere (|2.1ip . 

Discrete constraints. In order to determine the vector multiplier f we must specify 
discrete constraints on E(n) x E{n). A naive choice that imitates the form of the continuous 
constraints (12.91) is 



Ar k + pSl fc 7 = 0, 

where U k = R k+ iR\" G SO(n) is the finite rotation of the sphere in the space frame. However, 
this choice allows that the displacement Ar k of the center C may be not "horizontal" (i.e., 
orthogonal to the vertical vector 7), which is incompatible with the mechanical setting of the 
problem. 

For this reason, our choice of the corresponding discrete constraints (|3.4p will be 

Ar k + ^(n k - nl)i = 0, (4.9) 
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which does ensure that Ar k is orthogonal to 7. In the continuous limit (l k 1— ► I + u it transforms 
to the constraint (|4.2|) . 

Then, in view of (|4.9p and the second equation in (|4.7p . 

/* = -^(fl* - nftl + ^(^-1 - n T k _M 

and, therefore, 

F k Ml = -^[n k r k + r k n k ] + ^[njr fe + r fc n£] + ^[o fc _ir fc + r fc n fc _i] 
— ^-[^fc-i r fe + r fc^fc-i]> 

where, as in Section 2, T k = j k <8> 7^ = 7fc7fc 1 - 

As a result, we eliminate the multipliers -F fc and from (|4.7p . (|4.8p get the equations 

2 2 

2 2 
+ ^PVii* + r fe fi fc _i] - ^-[n^l^fc + r^j, (4.io) 

7fc = ^LiTfe-i, (4.11) 
which define an implicit map C: (£l k -i, 7it-i) 1— > (^fc>7fc)- 
Proposition 1. XTie map C admits the following compact representation 

K-k = ^fc-i^fe-i^fc-i, Tfc = M k _ 1 T k _iQ k _ 1 , (4.12) 

where JC k is the discrete analog of the momentum with respect to the contact point P of the 
sphere, 

K. k = n k (j + ^r fe j - (j + + ^(r k n k - n£r fc ) 

= M k + j(n k T k -T k nl) + ^(T k n k -nlr k ), D = m P 2 . (4.13) 

Indeed, in the continuous limit £l k i— > I + eu;, £<1 the matrix /C/e tends precisely to the 
angular momentum (I2.13P of the n-dimensional sphere with respect to its contact point, whereas 
the relations (|4.12p transform to the continuous equations (|2.12p . 

The map C therefore is evaluated as follows: 

1) given r^-i, determine lC k -\ by (|4.13p : 

2) Calculate JC k , T k by <KT2\i : 

3) Given IC k , T k , reconstruct Q k as a solution of the matrix equation (|4.13|) . 

The latter equation, in general, has several solutions, hence the map C is multi-valued. 
An an immediate corollary of Proposition [IJ we obtain 

Proposition 2. Regardless to a branch of the map C, it preserves the set of momenta integ- 
rals (|2.14p of the continuous Chaplygin sphere problem with K replaced by K,. 
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This property gives a solid justification of our choice of the discrete constraint (|4.9p . Note 
however that a discrete analog of the energy integral (|2.15|) in the generic case is unknown (see 
also Proposition [3] below). 

Proof of Proposition [TJ The second equation in (|4.12p follows directly from (|4.1ip . Now rep- 
lacing Tfc in the right hand side of (|4.10p by nk_ i r re _if2 re _i and using the identity fi re _ififc-i = I, 
we get 



mp 



2 L ^ ^ — re — 1 *v -l re — ij 2 

which, in view of (|4.13p . is equivalent to the first equation in (|4.12p . 



Remark 3. The form of equations (|4.12p . (|4.13j) is reminiscent to that of the discrete Euler top 
on SO(n) first described in |19j . 



M k = Ci^-i^-i. M k = n k j - ml, 



(4.14) 



however, the discrete momentum /C in (|4.13p depends not only on £1, but also on 7. Note that, 
in view of Proposition [TJ the subsequent momenta K, k -\,K, k admit curious intertwining relations 

jc k _t = Q fe _! U + j (r fe _! + v k ) \ -(j+ j (r fe _x + v k )\ 
Kk = ( J + § (IVl + " n *- 1 ( J + T (IVl + 

which are reminiscent to the following relation between the subsequent momenta of the discrete 
Euler top (|4~T4j) (see [16]) 

M fc _! = £l k -\J - ml_ l7 M k = Jf] fc _! - J. 

5 Discrete Chaplygin sphere on -E(3) and its particular solutions 

In the case n = 3 we can make use of the following parameterization of the body finite rotations 
/<?o + <?i - <?1 - <?| 2 (<7i<?2 - qm) 2(qm + 929o) \ 



2(9192 + 9390) ql + ql - q\ - ql 2(9293-9091) 

V 2(gig 3 -g 2 go) -2(9293 - 9o9i) 90 + ^3-91-92/ 



(5.1) 



where (/o> 9i are the Euler parameters subject to relations 9o + 9i + 92 + 93 = 1 an d 9o > 0. The 
operator (|5.1|) describes a finite rotation in M 3 about the vector q = (9i,92,93) T by the angle 6 
such that qo = cos (9/2) and |q| = sin(#/2) (see, e.g., [20]). 

Let now M k ,K, k £ M 3 be the vector representations of M k ,fC k G so* (3), 

M = (Mi, M 2 , M 3 ) T = (M 32 ,M 13 ,M 21 ) T , 

K, = (/Ci,/C 2 ,/C 3 ) T = (/C32,A^13,^2l) T - 

Then, in view of (I4.13p . we get simple expressions 
((■h + J?,)qoqi + {J2 - ^3)9293 \ 



M 



{J\ + J?,)qm2 + (J3 - Ji)qm 

\(Ji + J 2 )9o93 + (Jl ~ ^2)9192/ 



/C = M + 2Dg (q-(7,q)7)- 



(5.2) 
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In the continuous limit, when the angle 6 is small, 8 = euj, £ < 1, we have q 
qo = I — 0(s 2 ), and, up to terms of order e, 



w + 0(e 3 ), 



Ik 

£ 



Auj - D<7, w)7, A = diag( J 2 + J 3 + A + ^3 + D, J\ + J 2 + I?), 



which coincides with the expression for the angular momentum vector K in the classical Chap- 
lygin problem. 

As a result, in the vector form the map (|4.12p reads 



lC k = Q^^Kk-i, 7fc = njt_ 1 7fc_i, 



(5.3) 



where Qk is recovered from /Cfc, 7& by solving (15. 2p with respect to q and substituting the solution 
into (|5.ip . One can show that for real /C, 7, the equations (|5.2j) have at most 4 and at least 
2 real solutions q. 

Like in most of other discrete systems, in order to choose one of the branches of the map (|4,12p 
one should restrict to the case of sufficiently small finite rotations O^. In this case only one of 
the above real solutions q will be small and should be taken as the appropriate branch. 

Like the continuous system (12. 3p . apart from the geometric condition (7,7) = 1, the discrete 
system preserves two independent integrals 



(fC, 1 ) = h, {K,K) 



n. 



However, as simple numerical tests show, the energy integral (|2.5p is not preserved. 

The special case K, || 7. Like the continuous Chaplygin system, the map (|5.3p has the 
special case when fCk = h"fk, h =const. In view of (|5.2p . this implies the following relation 
between 7 and q 



[h + Dqoij,^, 



J 



^ 02 + Js)qoqi + {J2 - ^3)^2^3^ 

2 (Ji + J 3 )?o92 + (J3 - J\)qiqz 
\(J\ + J2)qmz + (Ji - h)q\q2/ 

Then one obtains the reduced map Gh '■ S 2 1— ► S 2 , Ghi'Jk-i) = 7fc evaluated as follows: 

1) given 7jt_i, one recovers q as a solution of (15. 4p . 

2) one calculates by (|5.1|) and 7^ by the second equation in (|5.3|) . 

Proposition 3. For any branch of the map Gh, it has the quadratic integral 
(7,A~ 1 7) =1. 

Proof. In view of (|5.4p . (|5.ip . and the second equation in (|5.3jl Fl. 

/ {Ji + Jz)qoqi + {J2 - J3)q2qs\ 



(5.4) 



(5.5) 



Ik 



2 



7fc+i 



{J\ + </ 3 )<?o<?2 + (^3 - J\)qiqz 
\{J\ + J2)qoqz + - J2)qiq2j 

( ( J 2 + J 3 )<?o<?i - (^2 
(^1 + </ 3 )<7o<72 - (^3 
\( J\ + J 2 )<? <?3 - (^1 



cr fc = [/i + £>(7fc,goq)]) 



^3)5253^ 

<^l)<7l<73 

J2)q\q2j 



(5.6) 



Substituting these expressions into (|5.5p . we obtain (7^., A 7^) = (7^+1, A 1 7/ c+1 ). 



3 To simplify notation, here and below we omit the discrete time index fc in the components of q. 
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As a result, in our special case the complex invariant variety of the map Qh is the spatial 
elliptic curve £, the intersection of the unit sphere (7,7) = 1 with the quadric given by (|5.5|) . 
The curve £ is 4-fold unramified covering of the planar elliptic curve 

£ = {w 2 = -(z - A^)(z - A^)(z - A^)(z - I)} 

and the points of £ are parameterized by the Jacobi elliptic functions associated to £q. Assume 
for concreteness that A^ 1 < A^ 1 < I < AT 1 . Then the parameterization reads (see e.g., [201 113]) 

71 = Ci cn(u|&), 72 = C 2 sn(u\k), 73 = C 3 dn(u\k), (5.7) 



where u is a complex phase parameter and k is the modulus of £q given by 

2 _ (A 1 -A 2 )(1-A 3 Q 
(A 2 -A 3 )(A 1 /-1)' 

Therefore, for a fixed I, the map Qh is reduced to one-dimensional map 

u k+ i = u k + Au k (u k , I), (5.9) 
with the increment Au k is a function of u k , I to be determined below. 

Remark 4. As seen from (|5.6p . the structure of the map Qh is similar to that of the Veselov- 
Moser discretization (j4.14j) of the Euler top on so(3). As shown in [161 [T]. t ne latter discretization 
preserves the momentum and the energy integrals, whereas the corresponding increment Au k 
does not depend on the argument u k in the elliptic parameterization. The same holds for another 
exact discretization of the top found in [7]. 

In comparison with the Veselov-Moser discretization, our map Qh involves the extra factor a k , 
which is not constant on the orbits of Qh- This gives an indication that under our map Qh 
the increment Au k depends essentially on u k . Moreover, as Qh is multi- valued, the function 
Au k (u k ,l) is expected to be multi-valued as well. This stays in contrast with the solutions of 
the continuous equations (|2.7p obtained from the Chaplygin sphere system under the condition 
K || 7: the components of 7 are elliptic functions whose argument changes uniformly with time t. 

In this connection the following natural problem of the complete integrability of the map Qh 
arises: is there a (multi- valued) reparameterization u = f(s,l), s £ C such that 

Uk = f(sk), u k+i = f(s k + As) 

and the increment As does not depend on s k ? Note that this property is equivalent to that the 
one-dimensional map (|5.9p preserves the invariant measure fi(u)du, where the density fi(u) must 
satisfy the discrete Liouville equation 

du k+ i _ n{u k ) 
du k n(u k+ i) 

and is related to f(s) by the formula df/ds = 1/fi. Due to the existence of the integral (|5.5p . 
this is also equivalent to that the map Qh itself preserves an invariant measure on S 2 . 

Since u is the argument of an elliptic function, the restriction of the map (|5.9p onto the 
real axis describes a diffeomorphism of a circle with the angular function Au k (u k , I). Hence, the 
answer to the above integrability question requires studying the properties of the diffeomorphism 
and applying known theorems in this field (see, e.g., [T7] and references therein). 

Determination of Au k (u k ,l). For convenience, in the sequel we omit the discrete time 
index k and denote Qh{l) = 7> u + A11 = u. To determine the function Au(u, I), we shall use 
the following property. 
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Proposition 4. The map Qh admits the following implicit form that involves the sum and the 

difference 0/7, 7 only: 

7-7 = j (7 + 7) x A- 1 (7 + 7 ), (5.10) 
2cr 

K = 



1- I (7 + 7, A- 2 (7 + 7)) 



h 

a = p . (5.11) 

1- ^<7 + 7,A- 1 (7 + 7)> 

Proof. As follows from (15.61). 



4(70 . 4 , . 

7 + 7=— Aq, 7 - 7 = -Aq x q (5.12) 
a a 

and, therefore, 

7 - 7 = Y2 (7 + 7) x A" 1 ^ + 7). 
Using (|5.12p . it is easy to check that 



2ql = 1 + A /l - ^(A-i( 7 + 7), A-!(7 + 7)) 



and that the factor <r defined in (15. 6p is the solution of the equation 

a = /i + Z^7,|A- 1 (7 + 7 )y 
Then, in view of relation 

(7,A~ 1 (7 + 7)) = i( 7 + 7,A- 1 (7 + 7 )), 

one obtains ([Q1j> . (I5TH . ■ 

Below we assume /i = 1 in (|5. 11 j) . which can always be made by appropriate rescaling of A 
and D. 

Implicit maps in the symmetric form (I5.10|) were studied in pQ as special discretizations 
of the Euler top on so(3). The symmetry allows to determine the function Au(u,l) in (|5.9p . 
Following [l], we set 

u+ u u — u 

u = , ou = . 

2 ' 2 

Then, due to addition formulas for the Jacobi elliptic functions (see e.g., [13]), the parameteri- 
zation (|5.7p . (|5.8p implies 

71 — 71 = -^-sn(-u)dn(-u)sn(5u)dn((5n), 71 + 71 = ^^-cn(n)cn(5u), 

X> = 1 - /c 2 sn 2 (u) sn 2 (5n) (5.13) 

and similar expressions for the other components of 7 — 7, 7 + 7. Substituting them into any of 
the equations (|5.10p . one obtains (see also [1]) 

2D = F cn(fa)dn(fa) y2 = (A 2 - A 3 )(AiZ - 1) 
x sn(5u) A1A2A3 
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or, in view of (|5.1ip . 

Hi = 4(7 + 7,A _1 (7 + 7)>, n 2 = ^(7 + 7,A~ 2 (7 + 7))- 
Then, making squares of both parts of (|5,14p and taking into account the expressions 

lfCf C$ C 2 2 \ 1/ 2 \ 1-A 3 / 



2?VAi Aa Aa^^V-pV ^ ^ A 2 -A 3 ' 

n 2 = i-(^cn 2 (n)cn 2 (5n) + %sn 2 (u)cn 2 (<5u)dn 2 (c^) + %dn 2 (u)dn 2 (5-u 



1 ( G _ v _ v ^ 2 {5u)An 2 {5u) 



V 2 \sn 2 (5u) sn 2 (5u) 



where 



G = -^cn 2 ^) - cn 2 (<5u)dn 2 (H + ^|dn 2 (,5n) = V 2 (l + asn 2 {5u) - (3sn 4 (5u)), 
k A 1 A 2 A 3 

= 2A!A 3 - (Ai + A 3 - A 2 ) AxA 3 1 - A 3 Z 

a (A 2 - A 3 )(AiZ - 1) ' P (A 2 -A3) 2 A!Z-1' 

after simplifications we obtain the following relation between the Jacobi elliptic functions of 5u 
and the function sn(tZ) (contained in T> only) 

—4- v2 (GV - V 2 cn 2 (5u)dn 2 (5u)) 

+ d l ~^ S ^ 6u \ (GV-V 2 cn 2 (5u)dn 2 {5u))-Vw{5u)V)+V{2Vw{5u)-G) = 0, (5.15) 
w(du) = sn(<5u)cn(<5«)dn(<ki). 

Note that in the case d = 0, when the map Qh reduces to the discretization of the Euler top, 
the terms with T> in (|5,15p can be eliminated and the increment Su becomes independent of u, 
as expected and was shown explicitly in pQ. 

Now, expressing sn(-u) =sn(u + Su) in terms of elliptic functions of u and Su, and using 
relations cn(u) 2 = 1 — sn 2 (u), dn(ti) 2 = 1 — k 2 sn 2 (u), we can rewrite (|5.13p as 

V = h U2 2,\ 2,* , 12 (fe 2 sn 2 (n)sn 2 (fa)(2 - (k 2 + l)sn 2 (Su) - k 2 S n\Su)) 
[1 — £rsn z (u) sn z (du)\ A 

+ 2k 2 sn 2 (Su)w(Su)W(u) + k 2 sn 4 (Su)[l - £; 2 sn 2 (u)] - 1), 
W(u) = sn(it)cn(it)dn(it) , w(5u) = sn(<5u)cn(<5u)dn((5u). (5.16) 

Substituting this expression into (|5.15p . one obtains a rather complicated equation for sn 2 (tt), 
sn 2 (5«), w(6u) that implicitly describes the dependence the function Auk{uk,l)- 

One can make this equation algebraic by introducing the variables p = p(Su), q = p'(Su), 
which are coordinates on the curve £ represented in the canonical Weierstrass form 

g 2 = 4(p- ei )(p-e 2 )(p-e 3 ), ei = ^(2-/c 2 ), e 2 = I(2A; 2 - 1), e 3 = -\{\ + k 2 ). 
Under the rational parameterization 

sn (du) = , cn (du) = , dn (du) = , w{d u ) = 

p - e 3 p - e 3 p - e 3 2 (p - ez) 6 
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equation (|5.15p together with (|5. 16j) defines a plane algebraic curve whose intersections P^' 
with £q give the values of 5u via the integral 



5u {j) 



dp 



2\7(p-ei)(p-e 2 )(p-e 3 ) 



6 Conclusive remarks 

In given paper we constructed a discretization of the Chaplygin sphere problem on E(n) which 
preserves all the momenta integrals and, in the particular case of the motion, the energy. This 
case is reduced to the multi- valued map Qh on so* (3) which is reminiscent to the Veselov-Moser 
discretization of the Euler top and is given by an addition law on elliptic curves. However, 
in contrast to known integrable maps of this type, the increment depends not only on the 
energy constant, but also on the point on the curve. According to [I], this implies that Qh does 
not preserve the standard Lie-Poisson structure on so* (3), which is quite expected due to the 
nonholonomic origin of the map. 

The integrability of the map Qh (which is interpreted as the possibility of writing its explicit 
discrete time solution) was shown to be equivalent to preservation of an invariant measure on 
the Lie algebra. Note that the latter property is often related to integrability in quadratures in 
the continuous case (recall the Jacobi last multiplier theorem). 

In this connection it would be interesting to develop integrability criteria for discrete systems 
preserving an invariant measure, in particular to formulate a discrete version of the Jacobi last 
multiplier theorem. 

Although the energy (|2.5p of the Chaplygin sphere expressed in terms of the momentum /C 
is not preserved in the discrete setting in general, this does not exclude the existence of another 
integral that transforms to the energy in the continuous limit. 

Lastly, as shown in [18], the classical Chaplygin sphere problem is a natural example of 
a nonholonomic LR system on the direct product SO (3) x M 3 (the kinetic energy is left-invariant, 
while the constraint distribution is right-invariant with respect to the group action). One can 
show that the discrete Lagrangian (|4.1I) and the discrete constraint subvariety given by (14. 9\\ 
possess the same properties, i.e., our discretization of the Chaplygin sphere is a discrete LR 
system on SO(n) x M. n . 

We then believe that it is worth studying properties of generic discrete LR systems and of 
their reductions, such as preservation of momenta, energy, and an invariant measure. 
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